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Hamilton— Jacobi's equation and Arnold's diffusion near invari- 
ant tori in a priori unstable isochronous systems. 

G. Gallavotti 
Fisica, Romal 

Abstract: Local integrability of hyperbolic oscillators is discussed to provide an intro- 
ductory example of the Arnold's diffusion phenomenon in a forced pendulum. 
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1. Introduction. 



Here I expose, in a simple case, a well established technique. I shall show that a 
Hamiltonian system composed of two clocks (i.e. isochronous rotators) and a hyperbolic 
oscillator can be integrated by quadratures in a full neighborhood of the unperturbed 
equilibrium of the oscillator. The rotators angular velocities form a vector lu — (u>i,u>2) 
' which we suppose Diophantine, see §2. 

The system is described by i — 1 pairs (A, a) £ R l ~ l x T £_1 of "actions" A = 
■ (A\,...,Ai) and "angles" (ax, ■ • • , a^_x) and one pair of conjugate coordinates (p,q) 

with a Hamiltonian: 



w • A + J(pq) + ef(a,p,q) (1.1) 

with J(x) and f(a,p,q) analytic for \x\ < k, \p\, \q\ < k 1 / 2 and a S T* -1 . We suppose 
that J'(0) > and that k > is small enough so that for some g > it is g Q < 
|Re J'(x)\ < \ J'(x)\ < 2g for all |a;| < k. 
' The Hamilton Jacobi's equation for the integration of the system is: 



u ■ d z $ + J((p' + d q S>)q) +ef(a,p' + q) - J(p'(q + d p <$)) (1.2) 



where the unknowns are $(a ,p', q), J(x), and we require that J — J be divisible by s. 
If $, J exist and are analytic in their arguments and in e, it follows that for e small 
enough there is a canonical change of coordinates transforming the Hamiltonian (1.1) 
into: 



u-A' + Jip'q') (1.3) 

hence the points with p 1 = q' = form a family of invariant tori parameterized by A'. 



And the motions are completely known, and very simple, in the new coordinates. 

a 

Theorem 1: the Hamilton Jacobi's equation (1.1) admits an analytic solution for alls 
small enough. The solution is analytic in e and in p' ', q, divisible by e and so is J — J . 

This (form of a well known) theorem is very useful because it provides full control of 
motions dwelling a very long time around the invariant tori, in spite of the nonlinearity 
of the equations. This turns out to be the key feature needed in the theory of Arnold's 
diffusion. 

In §2 I prove the theorem (following [CG] where, however, it is discussed in a much more 
general case). In §3 I introduce the notion of Arnold's diffusion, [A], in the example 
of a quasi-periodically forced pendulum (i.e. in a system consisting of a pendulum in 
interaction with two clocks) and discuss how to apply the theorem to show its existence. 
In §4, §5 I provide the details, from [G3]. 

The example may be useful as an introduction and to clarify the intrinsic simplicity of 
a problem that was solved by Arnold many years ago and that, outside a small circle of 
specialists, still seems to be (unreasonably) presented as a very hard problem. 
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§2. A classical perturbative analysis. 

Consider the Hamiltonian: 



Ha = UL ■ A + J (x ) + fo(a,po,q ) (2.1) 

with Jo holomorphic in Q Ko = {|xo| < «o} and /o holomorphic in P^ . Ko = {e~^° < 
| e laj | < , for j = 1 ,...£ — 1 } x { |po | , |<Zo \ < K o^ 2 } ■ Here £ is the total number of degrees 
of freedom. We assume that u enjoys a Diophantine property: \uj ■ v_\ > Co\v\~ T f° r 
some Co, t > and for all non zero integer component vectors v_ = . . . , vi-i) G Z 1 . 

def 

We call eo — ||/o||{ ,«o the lowest upper bound of |/o| in its domain of definition and we 
call go — HJqIUo the lowest upper bound to the derivative of Jo m its domain; we also 
need the greatest lower bound to Re Jo that we call g . We set T = min{g , Co} and by 
the assumptions following (1.1) it is g < |Re J'(x)\ < |J'(a;)| < 2g . Summarizing the 
constants eo,g ,T are defined so that: 

so = \\fo\k , Ko , 9o < \ReJ'(x)\ < \J'(x)\ < 2g , T = mm{g ,C } (2.2) 

The system has, therefore, two time scales, namely the pendulum time scale g _1 and 
the clocks time scale Cq 1 : T is the smallest. 

Let us define the "generating function" A 1 -a + piqa + < f ) o(a ,Pi, <Zo) via a function 
®o(a,p,q) = -(go(x) + (J - 9) _1 /o(«,P, q) where: 

def 

• we set go(x) equal to the derivative J' {x) of Jq{x): go(x) = J'q{x). 

• we set d = qd q — pd p and d d = d a . 

•fo(a,p,q) = h{a,p,q) - f (pq) and f (pq) is obtained from f(a,p,q) by taking 
the average over a and then expanding in powers of p, q the result and keeping only 
the terms that depend on the product pq. This means that we write f(a,p,q) = 

EZn=ofmAa)p m q n = Y, m , n ,v f ' « ,m ^ ^ P™ q" and Set J(x) = EJo,™,^. ^ 

this definition we keep control of the size of / because \f m ,n{(x)\ < £o^q 2(" l +") an j 
\f„, m ,n\ < e Ko l(m+n) e- €o|£l , where \ u\ = £\ N- 

Under the conditions on the auxiliary parameter 5q > specified below and in the 
smaller domain P^ _ SoiKoe -s , the generating map: 

po =pi +d go $ (a,pi,qo), qi = <7o + d Pl $ {a ,p 1: q ), A = A 1 + da®o(a,p 1 ,q ) 

(2.3) 

defines a canonical map Co : P^ s , Ko e- s « ~~ * P^a^a'- 

{Ao = Ai + Ki{a,Pi,qi) 
p =Pi+L 1 (a,p 1 ,q 1 ) (2.4) 
q = qi + Li(a,pi,qi) 

where H_ 11 Li 1 Li are suitable functions holomorphic in ^ _a K „ e -*o ■ 
In fact by simple dimensional estimates ("Cauchy estimates"), setting £ = £o — \ So, 
and k' = K e~? So and choosing B large enough it is (recall the above definition 

def 

£ o = ll/ollfo.Ko! see also appendix A2): 

\\d q $o\k' oK , ||0 P *o|le^, i* h \\d**o\\%,n' < BoTo^K-h^-^eo 
WWqM&kI < Bor - 1 Ko- 1 5o'" T " 2 eo 
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hence under the condition BiTq X Kq 1 eoS ~ i ~ T ~ 2 < 1, with B\ large enough, the map Co 
is defined (by the implicit functions theorem) and maps P^ _s 0tKQe -s into P^ , Ko - 
The choice of $o is just right so that by evaluating the Hamiltonian in the new coordi- 
nates it takes the form: 

Wi = uj ■ A l + Jiipm) + fi{a, Pl , qi ) (2.6) 

with Ji(x) = Jq(x) + f (x) and a suitable f\ "which is of order ||/o|| 2 "- 
In fact set £1 = £ -2<5o, «i = Koe~ 2S " &ndg 1 = g Q -B 2 e K^ 1 S^ 2 ,g 1 = 2g +B 2 SoK~ 1 6~ 2 , 
with B 2 large (see below). 

Note that B 2 eoKo 1 5 ~ 2 is a dimensional bound for 4-f(x) for \x\ < K e~ s °/ 2 (the S^ 2 
is due to the fact that the series over p, q is a double series: hence isolating the terms 
(pq) m and trying to bound the series expressing / (via the bounds on the Fourier-Taylor 
coefficients fv,m,n gives a factor (1 — e -5 "/ 2 ) -1 , and furthermore one has to differentiate 
/ with respect to x which yields a further division by k (1 — e~ s °/ 2 ). Then: 

9i < \ReJ{\ < \J{\ < 51, ex = \\fi\k uKl < B^^S, 2{l+T) ^ (2.7) 

1 o«o i 

provided B^EoTq 1 k^S^ 1 ^ 1 ^^ 1 < 1 for some B4 large enough and g 1 > ^g ,gi < 4c/ - The 
bound (2.7) can be checked by a first order expansion of the Hamiltonian Ho evaluated 
by expressing p 0} A l in terms of p\, qo, a via (2.3) and by developing to first order with 
respect to the increments po — pi, qi — qo, Ao ~ £i- The choice of $0 is what is needed 
to eliminate the first order terms and the remainders are easily estimated via Cauchy 
estimates, in a domain slightly smaller than the one in which Co is defined, which we 
recall to be P^ Q s Koe - s o ■ For details see Appendix A2 below (see also [Gl], Ch. V, §12). 

Proceeding as usual in the KAM proofs, [Gl], define £j = £ — 2(5 — ... — 2<5j_i and 
Kj = Koe- 260 --- 2 ^- 1 and fix 5, = £ ((j + 1)4)~ 2 . Then it is: k, > K e-^° and 
£j > |^o- We also set T 3 = min^-, Co) and g j+ i = g 3 + B 2 e j Kj 1 SJ 2 ^ +T ^ 4 . 

In this way under the conditions B^jT^ 1 k q 7 1 5J 1 ~ t ~ 1 < 1 and 2g > gj > gj > ^g a it 

will be: e j+ i < B^ 2 ^ 1 K^g^ 1 ^" 1 (j + 1)" having set q d = f 2(1 + r) + 4. 
Defining \j+i = X 2 (j + l)~ q for j > 1 and Ao = 1 we get Xj > e~ q ° 23 with a qo bounded 
proportionally to q. Furthermore defining, for a suitably large B 5 , r/j = B 5 g TQ 1 ^ q XjEj 
we find, from the definition of A 3 : 

Vj < Vo (2-8) 
Therefore, if B 6 is suitably large so that the condition: 

BeeoTo %To W < ^ q = 2(£ + r)+4 (2.9) 

encompasses all the ones found so far, one can define not only Co but also, recursively, 
C 3 in the domain P^, Kj , with image in P^ j _ l , Kj _ l . 

def 

The composition C = . . . o C 2 o C\ o C of the canonical maps will be well defined and it 
will give as a result a map C : ( A ,po, qo) i — ► (Aoo'-Poc loo) still of the form (2.4) with 
some function H_,L,L a 

into the "normal form" : 



some function H^,L,L analytic in _i 5o and transforming the Hamiltonian (2.1) 



Hoo = W ■ A^ + Joo(Poc<7oo) (2.10) 

with Joo holomorphic in Q _i (o wit < |Re J^J < [J^j < 47? . 
Comments: if fo(a ,p, q) = /(— a,q,p) then the symmetries of the problem imply that: 

H 00 {a,p,q)=H 00 (-a,q,p), L 00 (a,p,q) = L 00 (-a,q,p) (2.11) 
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a limit as j — ► oo of the corresponding relations which also hold identically for H ■ , Lj , , 
for j = 1,2,.... Furthermore has by construction average when evaluated at 
p = q = (a property of "twistless tori", see [G2]) or, more generally, at pq = 0. I shall 
call, therefore, A x the average actions of such torus. 

§3. An application. 

We consider a system with Hamiltonian: 

p 

n=u- A + — +g 2 (cos^- 1) +£/(¥>, a) (3.1) 

where (I, if) G R x T 1 describe a pendulum while (yl, a) G R 2 x T 2 describe two 
clocks and / is an interaction which we take to be a even trigonometric polynomial in 
a = ( Q!i, 0:2), ip. A non trivial case is / = cos(ai + ip) + cos(a2 + <p). 
Near the pendulum equilibrium one can use Jacobi's coordinates p, q, see Appendix Al. 
In such local coordinates the system has Hamiltonian: 

H'=w-A+J(jpq)+ef(a,p,q) (3.2) 

Hence by the theorem above we can find a (local) system of coordinates A ^ , a , poo , qoo 
in which the motion is described by a trivial Hamiltonian (2.10). This means that it is: 

A ca = const, a -> a + ut, Poc -> Pooe^ (x)t , gooe^ (:r) * (3.3) 

which shows that the set = q^ = is an invariant torus parameterized by A x and 
denoted T(A 00 ). 

Physically this is a very special set of motions in which the pendulum does not ever 
fall down from its unstable equilibrium position but performs small oscillations around 
it. It also shows that the sets q = 0, denoted W s (A OQ ), or p = 0, denoted W U (A 00 ), 
are invariant manifolds for t > and t < respectively, such that data on them evolve 
towards the (quasi periodic) motions on the invariant torus T(A OQ ), respectively, as 
t — > ±00. 

The meaning of w • A is the "energy of the springs that move the clocks" . TTie problem 
of Arnold's diffusion is whether one can find motions whose net effect is to transfer 
energy from one reservoir to the other for all e > no matter how small: note that if 
e = transfer is impossible. 

Mathematically this is a motion that starts near a quasi periodic motion with average 
action A = A and ends up, after a finite time, near one with average action A 1 7^ A . 
The key point in this definition is that Arnold's diffusion between A and A 1 if such 
motions exist for all e > small enough, no matter how small (provided £ ^ and, of 
course, for suitable e-dependent initial data). Here "near" means closer than the half 
the distance between the tori corresponding to i and A 1 . 

Of course for this to be possible it must be that the extreme motions have the same 
energy: it is a simple symmetry property of the above models (with / even) that the 
condition that two quasi periodic motions with parameters A 00 =A and A^ = A 1 
have the same energy is, generically and if e is small, simply uj_ ■ {A 1 — A ) = 0, see 
[G3]. 

Consider the line s — > A(s) = A + (A 1 — A )s, s G [0, 1], supposing it orthogonal to 
uj_ . One says that there is a heteroclinic chain between A and A 1 if there are M values 
So = 0, Si, . . . sm such that for all j's there exists a "heteroclinic" motion, i.e. a motion 
that is asymptotic as t — > —00 to the quasi periodic motion on the torus with average 
actions A x = A(sj) and that is asymptotic as t — > +00 to the quasi periodic motion 
on the torus with average actions A^ = A(sj+i). 

It is a simple application of the implicit functions theorem to show that any chain 
with Sj+i — sj small enough is a heteroclinic chain, see [G3], for generic /. In fact if 
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/ = cos(ai + p) + cos(a2 + f) this is the case, see [CG], [GGM]. The number Af will be 
in this case, and in general, of size O^ 1 ), see [GGM]. 

How small e has to be so that it is "small enough"? this depends on g. An interesting 
question is how small has e to be compared to g (the pendulum characteristic frequency). 

Consider a simple case: u = (n, 1) and g 2 = 77 and f{a,p) = [3 cos{a2+f)+er] c cos(ai + 
f). In this case if c is large enough taking (3 = ei] c is sufficient, see [GGM]. However one 
can do much better and allow [3 to be larger than 1, i.e. independent of e\ provided one 
excludes finitely many possibly exceptional values of (3, see [GGM]. This shows that a 
often stated condition that (3 has to be exponentially small in r^ x l 2 can be improved in 
three time scales problems. 

We now study the existence of diffusion reproducing for convenience of the reader §4 
and §5 of [G3]. 

§4. Geometric concepts. 

Let 2 k > be smaller than the radius of the disk in the (p, q) plane where the canonical 
coordinates ( A x , a ,Poo, <Zoo) can be used. I will drop soon the subscripts 00 to simplify 
the notation. To visualize the geometry of the problem we shall need the following 
geometric objects: 

def 

• (a) a point Xi, hctcroclinic between the torus parameterized by = A(si) = A i 
which we shall denote T(AA) and the "next" torus T(A i+1 ): this is a point Xi such 
that S t Xi f ^ +oc > T(A i+1 ) and S t Xi t ^_ oa > T{A i ). We denote the local coordinates of 
Xi as X, L = {Ai, ai,0,n). 

• (b) we can extend the manifolds W S (T(AA)) and W U (T(AA)) by applying to them the 
time evolution flow St- Then by our assumption the A i form a hctcroclinic chain and 
therefore W s (T(A i+1 )) evolves in a finite (negative) time to the vicinity of T(AA), and 
in fact its image (at a properly chosen negative time) will contain Xi. Hence it can be 
described in the local coordinates near such torus. The equations, at fixed q = k, of the 
connected part of W s (A i+1 ) containing X i} will be written in the local coordinates near 
^a,)Hs: 

Yi(a) = (A S i +1 (a),a,p s i+1 (a),K) (4.1) 

with I a — q_i\ < £ for some ( > (i-independent): it is A^ +1 (a i ) — A t , pf +1 (aj) = 
because we require ^(aj = Xi. There are constants F',F such that |il* i _ | _ 1 (a) — 
Ai +1 (a i )\ and max| a _ a ,\= fixed bf+i(s)| are bounded, for ( small enough, below by 
F'\a — a.i\ and above by F\a — a A; the constants F',F are generically positive, see 
[GGM] for instance. They are positive in the (non trivial) example / = cos(ai + ip) + 
cos(a 2 + ip). 

Note that W s (A i+1 ) also contains a part with local equations (A i+1 , a,p,0) which 
is not to be confused with the previous one described by the function Yi{a). This is 
more easily understood by looking at the meaning of the above objects in the original 
(A, a, I, ip) coordinates: in a way the first part of W s (A i+1 ) is close to p — and the 
second to <p = 2ir. They can be close to each other because of the periodicity of p, but 
they are conceptually quite different. 

• (c) a point Pi = Yi( a A) with | a ; - a J = r», where a^Ti will be determined recursively, 
and a neighborhood E>i\ 

Bi = {\A - A s i+1 {a)\ < Pil \a - aj < ft, \p s i+1 (a) - p\ < Pi , q = k} (4.2) 

where pi < Ti is another length to be determined recursively. If g d = ^<7 ,4c/ are a lower 
and upper bound to J'(x) for |a;| < 4k 2 , the point Pi evolves in a time Tj ~ g" 1 log« _1 
into a point X[ near T(A' i+1 ) which has local coordinates X[ = (A i+1 , a n, 0). 
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• (d) The points £ of the set Bi are mapped by the time evolution to points that, at the 
beginning at least, come close to T(A i+1 ) and in a time r(£) acquire local coordinates 
near T(A i+1 ) with p = k exactly: the time r(£) is of the order of g^ 1 logre -1 . 
If St is the time evolution flow for the system (1.1) we write S£ = <SV({)£ ( n °t e that S 
depends also on i). Then S maps the set Bi into a set SBi containing: 

Bl = {\A- A i+1 \ < 1^, \a - a'i\ < ^ Pi , p = n, \q\ < ^ Pl } (4.3) 

because all the points in B^ with A — A s i+1 (a), p = pf +1 (a), q = k evolve to points 
with A = A i+1 , p = k, q = and a close to a^, by the definitions. Here E is a bound 
on the jacobian matrix of S (which, being essentially a flow over a time Oilj^ 1 logK -1 ), 
has derivatives bounded i-independently: since we suppose that e is "small enough" we 
could take E = 1 + be for some b > if | A f - A l+1 \ < 0(e)). 

§5. The [CG]-method of proof of the theorem. 

Let TV be the number of elements of the heteroclinic chain. Consider the points 
Yf + i(a) £ W s (A i+2 ) with coordinates {Af +2 (a),a,pf +2 (a),K). They will evolve 
backwards in time so that A stays constant, a evolves quasi-periodically hence "rigidly", 
and pf +2 ( a ) evolves to n while the g-coordinate evolves from k to q = p s i+2 ( a ) (because 
pq stays constant, see (2.1)). This requires going backward in time by an amount of the 
order of T„ ~ g- 1 log^b t s +2 (a )T 1 + °o. 

Therefore there is a sequence a n ^ <2j+i such that a" — > a j+ b p|+2(ffi") ~~ * 
A 8 i+2 (a n ) —> A i+1 and a n — uj_Ta n n^oa ) as a consequence of the Diophantinc 

def 

properties of w. So that there is a i+1 = a n with n large enough and a point P%+i = 
(A S i+2(°L i+1 ), a i+1 ,p s i+2 (a i+1 ),n) £ W s (A i+2 ) (actually infinitely many) which evol- 
ves, backwards in time, from to a point of B[. 

Hence we can define r^+i = | a i+1 — a i+1 \ and pi+i small enough so that the backward 
motion of the points in Bi + \ enters in due time into B[. It follows that the set Bi evolves 
in time so that all the points of B i+ i are on trajectories of points of Bi. Hence all points 
of B_\f will be reached by points starting in Bq. 

This completes the proof. All constants can be computed explicitly, even though this 
is somewhat long and cumbersome. The result is an extremely large diffusion time T 
(namely ~g~ l times the value at W of a composition of J\f exponentials! at least this is 
the estimate I get after correcting an error in §8 of [CG]: the error was minor but its 
simple correction leads to substantially worse bounds, see [CG] ) . 

§6. Concluding remarks: 

(1) The above proof of Arnold's diffusion is very simple but it leads to a "bad" bound. 
One can improve the bound by using that there are many invariant tori and one can 
choose around which to construct a diffusing trajectory (for e small enough). The bound 
goes down to g- 1 0{2 u ), sec [G3]. 

(2) This is still far from the best bounds in the literature for similar problems, [Be], [Br], 
obtained by different methods. 

(3) The above (obvious) proof is likely to be what Arnold, [A], had in mind when he 
proposed his example: the example has in common with the above (1.1) the key feature 
of having a coordinate system like the one discussed in §2. The corresponding local 
Hamilton Jacobi's equation admits a solution "with no exceptions" , one says a "gap- less 
system of coordinates" , trivializing the flow: usually when the system is not isochronous 
or it is not of the Arnold's type there are open regions of phase space where the nice 
coordinates in which the flow is trivial cannot be defined, see [CG] for what can be done 
in such cases. 
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(4) The role of the special coordinates (A^, a,Poc, Qoc) defined after (2.9) is essential: it 
is clear from the analysis of §5 that one must control the motions near the tori over a very 
long time span. Since the equations of motion in the original coordinates are non linear 
this would be impossible in absence of an exact and simple solution, as the one given 
by (3.3): any small perturbation of the initial data would be amplified exponentially in 
time in an uncontrollable way. The (3.3), instead, confine the expansion to either the p 
or q coordinate and in a very explicit way. Hence the above is a good example of a non 
trivial use of the Hamilton-Jacobi's equation. 

(5) The method discussed in §5 is a simplified version of a method known as "windowing" 
described in a early work, [Ea], and developed in more recent works [M], [C]. But it is 
much less ambitious and detailed and the bounds obtained in the last two papers are far 
better than the ones we achieve here. I am indebted to P. Lochak and J. Cresson for 
pointing out the relevance of the latter papers. 

Acknowledgments: This is a contribution to the ISI school on dynamical systems 
of June 1997 in Torino: I am indebted to G. Gentile and V. Mastropietro for many 
discussions and suggestions. I thank A. Ambrosetti and E. Vesentini for the invitation 
to lecture at ISI. 

Appendix Al. Jacobi's map. 

This appendix is standard: here it is taken from A9 of [CG] with small changes, to use 
it for future references. 

The theory of Jacobian elliptic functions shows how to perform a complete calculation of 
the functions, below denoted R, S, in terms of which the canonical Jacobi's coordinates 
are defined, see [GR] (9. 198), (9. 153), (9.146), (9.128), (9.197). The result, reported for 
completeness, is discussed in terms of the pendulum energy: 

^-+ 5 2 (l-cos<^) = E (Al.l) 

where the origin in <p is set at the stable equilibrium, to adhere to the notations in the 
theory of elliptic functions. Setting u = t(E/2) 1 / 2 = e 1 ^ 2 gt, k 2 = 2g 2 /E = e _1 where e 
is the dimensionless energy so that e — 1 is the separatrix, let: 

r /2 da 

K ( fc ) = / n ^ ■ 2 u/2 ( AL2 ) 

Jo (1 — k sin a) 1 -'^ 

We shall use the "standard" notations (i.e. those in [GR]) for the Jacobian elliptic inte- 
grals except for x(.), which is usually denoted q(.), but which we would confuse with the 
variable q that we want to construct: 

_ 11 1,1/2 

k'-(l-k 2 V/ 2 a, -a - -e 1 / 2 a A = -- - 

fc (1 fc) , gj 9 2kK{k>) e 9, 21 + k 1 / 2 (A1.3) 

x(k') =e -*K(*)/K(*') = A + 2A 5 + 15A 9 + 150A 13 + 1707A 17 + . . . 

In terms of the above notations we have, directly from the definitions (i.e. from the 
equations of motion): 

I(t) = tp= -2ge 1/2 dn(u,k), ip(t) = 2 am (tge 1/2 ) (AlA) 

which yield, changing the origin for ip to the unstable point to conform with our notations 
(i.e. obtaining tp(t) = 2( am (tge 1 ' 2 ) + tt/2)), for I(t) = R(p(t),q(t)),ip(t) = S(p(t),q(t)): 

„ n , l0 dn (iu, k') . S 1 S .sn(iu.k') 

R=-2ge 1/2 ) ' , ' , sin- = — , cos - = i ), ' ( (A1.5) 

y ca(iu,k') 2 cn(m,fc') 2 cn(iu,k') ' 
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Setting p = e~ 9jt ,q = x(k')e 9jt , see [GR], and using R(p, q) = g,j(—pd p + qd q )S(p 1 q) to 
evaluate S from R, the quoted basic relations between elliptic integrals imply immediately 
that the I(t) = R(p(t),q(t)),tp(t) = S(p(t),q(t)), solve the pendulum equations if: 



R(p, q) = - 2 5J 



P i Q V^°° I 1 \n 1+z ( n 2 

1 (P 2 " 



S(p,q)=2 arctgp- arctgq - £~ =1 (-1)« i±|^ 



- 1 ) 



2n-l 



sin 



cos 



S(P,g) _£J 

2 9 
g(p,g) _ _ 



P 2 V^ 00 I i\n l-t" V „2rt-1 „2n-l> 



1-P 2 



1+p 2 



feS + 2ESLi(-i) 



n 1— a 2 " , 
1+rr 2 



(A1.6) 



with a; = pg. Note that </j depends on x, and so do k',k: hence the coefficients of 
the first and of the last two of (A1.6) are also functions of x = pq. Furthermore the 
(dimcnsionless) energy e becomes a function e(£) of £ — pq defined by inverting the map: 



£ = x(k') = jc((1 



-l ) l/2 ) 



(A1.7) 



and the point corresponding to ip = tt and to a dimcnsionless energy e, has coordinates: 



p=l, q = x(k') 



(A1.8) 



(a rearrangement of (A1.6) showing convergence for p = 1 and \x\ < 1 is exhibited below). 
The variables (p, q) defined above are nice and natural: however they are not canonically 
conjugated to (I, p): the Jacobian determinant of the map (p, q) < — > (I, ip) is not 1. But 
the Jacobian determinant must be a function D(x) = gjj'fy of x alone (i.e. of the product 

pq); then (Al.8) and the equations of motion imply that D(x)^ 1 — gj 1 2g = ^g de J x ■ 
Therefore one can modify the variables p, q into new variables {jpj, q,j) = (pF(x), qF(x)) 
with F such that the Jacobian d{ ^H ) = d( * { J p q q f D(x), which is D(x) ■ d x (xF 2 (x)), is 

^l/2^ e( a: ) 1 - /2 ~l ^l/2 



identically 1. One finds: F(x) = (4g) 
To invert the map (pj, qj) 

then: p — pjG(xj) and q = qjG(xj), x — xjG 2 (xj). The final result is a local canonical 
map between Jacobi's coordinates (p.j,q,j) and global (I,<p) coordinates: 



(pF(x),qF(x)) define xj d = pjqj and G(x.j) d = F(x) 1 



I = R(pjG(xj),qjG(xj)), p = S(pjG(xj),qjG(x.,)) 



(A1.9) 



where R, S are defined above, see (A1.6) which are written in a form easily recognized in 
the elliptic functions tables. The functions i?, S can be rewritten in the following form: 



x q 



1+x 2 



R(p, q)=- 4 9 [£^ =0 ( T - f |> 
S(p,q) =4 [Em=o( arct g x? >- arctgx m g)] 



sin 



cos 



S(p,g) _ 2g.,(x) 
2 9 

S(p,q) _ gj(x) 

2 2g 



x p 



x q 



Z^m=0V x ) \l+x 2m p 2 1+x 

1 — 9 Y^ 00 (— 1 \ m ( x2m p 2 _i_ x ""q' \ 

1 Z 2_/m=0^ 1 ) \l+x' 2 ™p 2 + l+x' 2m q 2 ) 



(ALIO) 



exhibiting some of the properties of the Jacobi's map in a better way. 

One checks that in the (pj, qj) variables the pendulum Hamiltonian, (Al.l) has become 
a function J(pjqj) = 2g 2 + gx j + 0(x 2 j) . The domain of definition of the map is given by 
the properties of the elliptic functions or, more restrictively, by the domain of convergence 
of the above scries. It includes a disk of some radius pj > around the origin. 

The important symmetry R(p,q) = R(q,p) and S(p,q) = —S(q,p) is manifest. 
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Appendix A2: Dimensional bounds. 

We give here two examples of dimensional estimates. 

• The bounds (2.5) are obtained by writing (with v e Z 1 ^ 1 ): 

p iv_- a m n 

*o(a,Pi,«b) = - E f"W- T / u 7 421 - 1 ) 

■ v +g(piqo){n - m) 

\v\ + \'n-m\>0 

so that setting £ ' = £ — X' K o = K - e ~~^ an d Co = Co — K o = K o e ~~^ (assuming I > 3 
and using the bound on \fv,m,n\ preceding (2.3)): 



l*o|| C » l(t »< 



E 

v ,m,n 



AUI^U + ^X_, 5o|ide _ i5o( „ +m) £Q < ^ £0^ 



for a suitably chosen B\; hence (2.5) follows because, for instance, |9 g $o I can be 

_ i 

bounded above by const k 2 ||$o||^',k^' (a typical example of a dimensional estimate). 
• By using (2.3) one immediately finds: 

Tin = Hi + (uj_ ■ 9$ + go(pm)0$o) + fo(a,p 1 ,q ) - f (pi qo)+ 

+ [M(pi + d qo $ )qo) - Jo(piqo) - Jo(piqo)qod qo $ ]+ ^ 

- [Jo(pi(qo-d Pl $o))-Jo(piqo) + Jo(piqo)Pid Pl ®o\ + 

+ {fo(a,pi + d qo <P ,qo) - fo(a,pi,qo)} - {J (pi (qo + d Pl $ )) ~7o(Pi *>)} 

The terms in square and curly brackets can be bounded dimensionally: for instance the 
first term in square brackets is bounded in the domain P^ -2S , Ko e- 25 o by a bound on the 
second derivative of Jo m P^ -s .Koe^ 5 a (*- e - 2 5o times Kq 1 (l — e _l5o ) _1 times a bound 
(on P^ _2s . Ko e- 2S o ) °f the square of pid Pl $>o or qod qo &o that is given by (2.5)': provided 
5q in P^ -28 ,K e- 25 o which gives the condition following (2.7)). 

Likewise the terms in curly brackets can be bounded. The final result is (2.7). 
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